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"ti ' Abstract. - In the context of the energy landscape description of supercooled liquids, we 

propose an explanation for the different behaviour of fragile and strong liquids. Above Gold- 
en , stein's temperature T^ , diffusion is interpreted as a motion in the phase space among saddles of 
the potential energy. Two mechanisms of diffusion then arise: mechanism A takes place when 
the system overcomes potential energy barriers along stable uphill directions, while mechanism 
■ O ■ B consists in finding unstable downhill directions out of a saddle. Depending on the mutual 
G ' efficiency of A and B, the usual classification of liquids in fragile and strong is recovered. 
O , Moreover, this scenario naturally predicts the possibility of a fragile-to-strong crossover when 
lowering the temperature. 



> 

■T^j- I After the seminal paper by Goldstein in 1969 M, it has become customary to think of the 

OQ . dynamical evolution of supercooled liquids in terms of motion of the state point of the system 
^— ^ I upon its rugged potential energy surface M. More precisely, at low temperatures, but above 
the glass transition, the diffusion of a system at equilibrium can be interpreted as the result 
Qs ' of two different processes: the thermal relaxation into basins defined by the many minima of 
the potential energy (mira-basin relaxation) and the hopping from basin to basin by crossing 
d ' potential energy barriers (mier-basin relaxation) ||l||. 

Crucial condition for this description to be correct is that the two relaxation times Tintra 

^_^ ^ and Tinter are well separated, that is Tinter S> Tintra- if these time scales are of the same 

Zh ■ order, it is not sensible to discriminate between the thermalization inside a minimum and 

Q ' the hopping among different minima. Indeed, the very requirement of the separation between 

O . these two time scales led Goldstein in Q to define and estimate a crossover temperature T^, 

above which this hopping-relaxation description is no longer valid, since Tintra and Tinter are 

not well separated. Below T^, on the other hand, crossing of potential energy barriers by 

/\ ' thermal activation becomes the primary mechanism of diffusion. The Goldstein temperature 

^ . Tx is in general higher than the glass transition temperature Tg-. in the interval Tg < T < T^ a 

■ - - ' supercooled liquid is still in equilibrium on experimental times and can be described in terms 

of relaxation into basins and thermally activated hopping among them. 

Some authors 0,0 have subsequently interpreted T^ as the temperature below which the 
ideal Mode Coupling Theory (MCT) [^ breaks down. This is because MCT is considered 
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unsuitable for describing activated dynamics, so that it prescribes a spurious dynamical tran- 
sition at the point where actually barriers crossing becomes a dominant process. According 
to [0,0] the temperature Tc where MCT locates this transition must therefore be identified 
with Tx ■ 

A useful realization of Goldstein's scenario has been introduced by Stillinger and Weber H : 
the phase space is partitioned into basins of attraction of the potential energy minima and each 
dynamical configuration is mapped onto the relative minimum, giving rise to a purely inter- 
basins dynamics. By comparing this pseudo-dynamics with the original one it is possible to 
directly verify the validity of Goldstein's decomposition into intra and mier-basins relaxation 
and estimate T,j. B. Computationally, the Stillinger- Weber scheme is realized by steepest- 
descent quenching of configurations generated by molecular dynamics simulation, and it proves 
remarkably useful in describing deeply supercooled liquids. 

The purpose of this Letter is to generalize Goldstein's ideas by introducing an energy 
landscape description of supercooled liquids which is valid also above T^, and, in so doing, to 
reach a better understanding of the physical difference between fragile and strong liquids ||3| . 
Within the Gibbs-DiMarzio entropic theory |8|, it is in fact possible to check a posteriori 
that fragile liquids must have a Kauzmann temperature Tk |9| close to the glass transition 
temperature Tg, while for strong liquids it must be Tr- '^ 0. However, this explanation of 
the difference between fragile and strong liquids has the drawback to involve the Kauzmann 
temperature Tj^^ whose existence is not universally accepted. Moreover, the fragile and strong 
phenomenology occurs at temperatures higher than Tr- and Tg, so that an interpretation based 
only on these two quantities seems unsatisfactory. 

At high temperatures Goldstein's description breaks down because there is no clear sep- 
aration between vibration inside a minimum and hopping among different minima. In order 
to find a new description we have therefore to start from this key observation: the hopping- 
relaxation scenario breaks down above Tx because the system is no longer spending most part 
of the time vibrating around minima. This happens because by raising the temperature, and 
thus the average potential energy, the system starts spending most part of the time in regions 
of the landscape where minima are rare. On the other hand, due to the great complexity 
of the potential energy landscape, other stationary points exist and are more numerous at 
such higher energy levels: these objects will not be minima in general, but saddles with some 
unstable directions. A generalization of Goldstein's description above T^ must therefore deal 
with saddles. Before giving a better specification of these objects, let us note that the rele- 
vance of saddles for the dynamics of glassy systems has been already recognized in the past 
in the context of mean-field spin-glasses |10|-[1^ . 

Saddles are stationary points of the potential energy with an arbitrary number K of 
negative eigenvalues of the Hessian, i.e. of unstable directions. The index K can take any 
integer value from K — Q (minima) to K = D (maxima), where D is the dimension of the 
phase space (for simple liquids D — Nd, where N is the number of particles and d the space 
dimension). The index density is defined as A: = K/D, and the potential energy density as 
u = U/N, where U is the extensive potential energy. It is useful to introduce the average index 
density k of the saddles located at a given potential energy level u. An important feature of 
this function k(u) is that it is in general an increasing function of u: the higher the energy, the 
larger (on average) the number of negative eigenvalues of a saddle at that particular energy. 
This fact has been explicitly proved for a Lennard- Jones liquid in |l3[ , where k{u) has been 
found to be a monotonously increasing function. 

When considering generic stationary points, it can be introduced a generalized notion of 
potential energy barrier, defined as the extensive energy difference AU between a saddle of 
index K and a higher, but adjacent, saddle of index K + 1. For K = this is the standard 
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definition of energy barrier between a minimum and a simple saddle. Provided that we know 
the average index as a function of the potential energy, fc(u), we can give a rough estimate of 
the average barrier A[/ at energy u by evaluating the energy difference between saddles with 
index difference IS.K = 1 (see also |TQ]). We have 

1 = AX-/c'(u)AC/d , (1) 

and therefore 

AU{u)^[k\u)d]-^ . (2) 

This estimate states that the potential energy barriers between stationary points at energy 
u change according to the change in slope of the key function k{u). Note that if there is an 
energy uq such that k'{ua) — 0, then barriers diverge at uq. This happens for example in the 
p-spin spherical spin-glass model: here both k{u) and k'{u) vanishes at the so-called threshold 
energy and a purely dynamical transition occurs (see, for example, |14| ). This is in agreement 
with the mean-field nature of that model, which requires the barriers among minima to be 
infinite. On the other hand, in the Lennard- Jones case k(u) has been found to be a linear 
function, always having nonzero slope |13|] : this implies that in such a system barriers among 
saddles are finite and independent of the energy level. 

A generalization of Goldstein's scenario may now be pursued by making the following 
hypothesis: the diffusion of a supercooled liquid can be described as a motion of the system in 
the phase space among the neighbourhoods of saddles of the potential energy. In order to give 
substance to this statement it is important to specify how to define the neighbourhood of a 
saddle. This can be done by considering the effective potential W = (VC/)^, i.e. the modulus 
square of the force: all the saddles of the original potential energy U are absolute minima of 
this effective potential W ||l^, |5| . At low temperatures the average potential energy of the 
system is small enough for the index k to be on average zero: the system stays mainly around 
minima and Goldstein's description is recovered. Indeed, below Tj. a quench of a dynamical 
configuration on the W surface will drain, on average, to a minimum of U, and the present 
scheme reduces to the Goldstein/Stillinger- Weber scenario ||l^. However, if the temperature 
is larger than Tx, a quench will typically lead to a saddle of U, meaning that Goldstein's 
scenario breaks down. 

Once assumed that diffusion can be described as a motion among the neighbourhoods of 
different saddles, we can identify two mechanisms for this motion to take place p7| : 

• Mechanism A • It consists in the crossing of potential energy barriers: the system, 
initially in the neighbourhood of a saddle, takes an uphill stable direction to reach a saddle 
with higher energy, eventually going downhill to a final saddle at roughly the same initial 
potential energy. Note that, even though we are considering global stationary points of the 
whole potential energy, a transition over a potential energy barrier will be local in real space, in 
that the rearrangement process will involve only a finite number of particles in the system |i| , 
with all the other particles acting as thermal bath. In this respect, mechanism A is compatible 
with both a canonical and a microcanonical description of the system. Note that potential 
energy is stored at the crossing point. At low enough temperatures mechanism A is driven 
by activation and its efficiency directly depends on the temperature through the Arrhenius 
transition probability 

Pa -^ exp{-AU/KT) , (3) 

where k is the Boltzmann constant and AC/ is the potential energy barrier, which can be 
estimated using the function k{u). Of course, the efficiency of mechanism A decreases by de- 
creasing the temperature. Were this the only mechanism of diffusion (as in a one-dimensional 
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phase space) , a knowledge of the exact form of the barriers as a function of the energy would 
be sufficient to predict the behaviour of the system. 

• Mechanism B • It exploits the fact that saddles have in general a non-zero index 
K: in this case the system can find an unstable downhill direction, which brings it out of the 
basin of the initial saddle. No activation is needed for this to happen. The system arrives at 
a lower potential energy level with an excess of kinetic energy, which is expended by climbing 
up again to a new saddle at the same initial energy. This mechanism takes place when a 
given cluster of particles suddenly finds a local rearrangement (an unstable direction), which 
sharply decreases their potential energy without any need of crossing a barrier. The extra 
kinetic energy acquired by the region is eventually dissipated by interacting with the rest of 
the system. Note that, unlike for mechanism A, in this case kinetic energy is stored at the 
crossing point. Mechanism B is the true signature of the multidimensionality of the phase 
space (as opposed to the typical one-dimensional picture of barriers hopping) and its nature 
is entropic, not energetic, with no direct dependence on the temperature. However, there 
is a crucial indirect dependence of mechanism B on the temperature, due to the fact that 
the average index k (ruling the efficiency of this mechanism) is an increasing function of the 
potential energy of the system, which is in turn an increasing function of the temperature T. 
For example, at very low temperatures the energy of the system is so small that only minima 
are visited on average, and of course mechanism B is frozen. 

Summarizing, the efficiency of both mechanisms decreases when the temperature is de- 
creased. A reasonable assumption is that, at a given temperature T , diffusion is ruled by the 
most efficient of these two competing mechanisms. Therefore we must compare the efficiency 
of A and _B as a function of the temperature, in order to understand which one of them drives 
the slowing down of the system. 

As a first step in this direction we introduce the threshold energy Uj/j, defined by the 
following relation: 

k{uth) EE . (4) 

Below uth unstable saddles become very rare and minima dominate, such that mechanism B 
cannot work. Note that, according to relation (0), energy barriers between threshold minima 
can be estimated from the slope of k{u) at the threshold. As already remarked, this slope 
is nonzero in a Lennard- Jones system p^ , giving finite barriers among threshold minima. 
Furthermore, in the Lennard-Jones case it has been observed that the threshold energy lies 
well above the energy of the lowest glassy minima found in the system [|3[ . 

A key feature of the threshold energy utt is that it allows us to define a critical tem- 
perature for mechanism B. To this aim we have to note that uth is the energy density of 
the threshold minima without considering the vibrational contribution. The total equilibrium 
average potential energy density of a system vibrating around a generic minimum with energy 

Urn, can be estimated as 

3 

u{T) ^Um + -kT , (5) 

where 3/2kT is the vibrational contribution within the harmonic approximation. In this way 
we can define a critical temperature Tb for mechanism B via the following relation: 

3 

UeqiTs) - -zHTb = Uth , (6) 

where UeqiT) is the equilibrium average potential energy density of the system. Thus, Tb is 
the temperature below which the system vibrates predominantly around minima, rather than 
saddles, and mechanism B is frozen. In other words, above Tb the typical saddle has got 
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an extensive number of unstable directions, that is K ~ 0{N), and mechanism B is highly 
efficient. At Tb an entropic bottleneck is created, for K is no longer extensive and the time 
needed to find an escape direction from a saddle diverges. It is important to understand the 
following point: at low temperatures mechanism A may be very slow, but in principle it is 
always available thanks to thermal activation. On the contrary, there is no activated regime 
for mechanism B: when the typical value of k is zero, the system cannot borrow an extra 
direction to escape a minimum. For this reason, mechanism B must freeze much more sharply 
than A, passing from an efficient phase {T > Tb), to a frozen one (T < Tb), with no activated 
intermezzo. 

In order to compare the two mechanisms we have to ask: What is the efficiency of mecha- 
nism A when mechanism B dies out, that is at Tb ? To answer this question we must compare 
the size of the potential energy barriers at Tb to the amount of thermal energy available to 
activation at this temperature, i.e. kTb- The barriers at Tb are given by the potential energy 
difference between threshold minima and simple saddles, and their size can be estimated via 
equation (||) as 

AU{uth)-^[k'{uth)dr' . (7) 

As already noted, in a Lennard- Jones system AU{uth) is finite p3| . We have no a priori way 
to know which one of the two quantities, kTb or AU(uth), will be the largest one. For the 
sake of simplicity we will consider the two extremes cases where one quantity is much larger 
than the other, classifying liquids into two groups: 

• Class I - The first class of systems is defined by the relation 

AU{uth) < kTb . (8) 

The potential energy barriers at Tb are very small, relatively to the thermal energy, so that 
mechanism A is still very efficient at the temperature where mechanism B freezes: activation 
is not even needed to overcome energy barriers at Tb, for too large is the kinetic energy of 
the particles. When the temperature is decreased further below Tb mechanism A (the only 
still available) slows down, and eventually Goldstein's temperature T^ < Tb is reached, where 
At/ ~ kT: in order to pass from minimum to minimum thermal activation is now needed 
to cross potential energy barriers, whereas B is completely unavailable. In the simplest case 
where potential energy barriers do not strongly depend on w, we expect the relaxation time 
to increase according to the Arrhenius law, until the glass transition Tg is eventually reached 
at lower temperatures. For Class I systems thus, the slowing down is entirely driven by the 
slowing down of mechanism A, i.e. by thermal activation. 

• Class II - For the second class we have 

AUiuth) > kTb . (9) 

In this case a very different behaviour may be expected: potential energy barriers are very large 
at Tb and therefore mechanism A is already very slow at the temperature where B becomes 
unavailable. Above Tb the system has no need to overcome potential energy barriers, because 
it can use a faster, non-thermally activated mechanism of diffusion, that is B. However, 
by decreasing the temperature Tb is eventually reached and mechanism B suddenly freezes. 
Thus, when at Tb for the first time potential energy barriers must be overcome, because no 
other mechanism of diffusion is available, these are already very large compared to kTb- In 
this case, therefore, there must be a sharp increase in the relaxation time, entirely driven by 
the slowing down of mechanism B. Note that in this case T^ = Tb, because only below Tb 
activated barriers crossing remains as the only mechanism of diffusion. Finally, if barriers are 
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indeed very large at Tb , the glass transition must occur close to Tb and activation must play 
little role in it. 

The smooth, uniform slowing down according to the Arrhenius law described for class 
I corresponds to the strong liquid behaviour, while class II exhibits a sharp dependence of 
the relaxation time on the temperature which is typical of fragile systems. Note that no 
energy scale, or other dimensional parameter, is involved in mechanism B, which has a purely 
geometric nature. This fact suggests that the increase of the relaxation time caused by the 
slowing down of B in the fragile case should be well reproduced by a power law. This point 
certainly deserves further investigation. 

The validity of the proposed scenario can by directly tested through equations (0), (Eh, (ph, 
(P) and (P), provided that the the average index k{u) as a function of the energy is known. 
This function can be easily computed numerically in models systems of liquids by sampling 
saddles of the potential energy. This program has been explicitly carried out in |13| ] for a 
Lennard-Jones liquid, where it has been found 

AU{uth) -^ 10 kTb . (10) 

According to the classification given above, relation (pl]| ) is typical of Class Il-fragile systems. 
This conclusion is in agreement with the fragile nature of Lennard-Jones liquids. Furthermore, 
in PS I it has been proved that Tb coincides with the MCT critical temperature Tc- According 
to [H3] this implies Tg ~ T^,, in further agreement with our description for Class II systems. 
These facts strongly support the validity of the present description, at least for the fragile 
case. The same kind of investigation should be performed in different systems, to see whether 
equation (^ is indeed a key feature of fragile liquids. 

There is an important phenomenon, whose existence can be predicted as a simple corol- 
lary of our scenario. Consider Class II, but suppose that barriers are not very large at Tb- 
According to what stated above, at Tb the mechanism of diffusion switches from B to A, and 
if barriers are now not too large here, mechanism A will not be completely unavailable yet. 
For example, we can imagine that the viscosity prescribed by activation at Tb is 

77a(Tb) ^ exp[AC/(ut,,)/AcTB] ~ 10**P , (11) 

quite far from the value rj ~ lO^^P defining the glass transition. What happens is then the 
following: Approaching Tb there is a very sharp jump of r/ (driven by the slowing down of B) 
up to the value r/ ^ lO^P. At Tb mechanism B is outstaged by the more efficient A and the 
viscosity starts increasing in a strong, Arrhenius-like way according to mechanism A, until 
eventually the glass transition occurs at lower temperatures, when rj ^ 10^'^P. A system of 
this kind must therefore display a fragile high temperature phase for T > Tb and a strong 
(Arrhenius) low temperature phase for Tg < T < Tb, that is a fragile-to- strong crossover. 
Note that the opposite strong-to-fragile transition cannot happen: if barriers are small at Tb 
the system cannot switch back to mechanism B at lower temperatures, because there is no 
activated regime for B. 

Remarkably, there are some recent evidences that a fragile-to-strong crossover as the one 
described above is actually present in some systems. Preliminary experimental observations 
indicate that this phenomenon occurs in supercooled water [Q, while it has been quite un- 
doubtedly observed a fragile-to-strong crossover in viscous silica, both from extrapolations 
of experimental data M] and from molecular dynamics simulations pO| . These are further 
positive tests for the validity of the scenario proposed in this Letter. 

To conclude, we briefly compare our description with the Instantaneous Normal Mode 
(INM) analysis of supercooled liquids pi[. A fruitful idea in the INM context is to relate 
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diffusion to the average fraction of negative eigenvalues of the potential energy Hessian |23| : 
instantaneous unstable modes are interpreted as the contribution of regions with negative 
curvature sampled by the system in crossing potential energy barriers. In our description 
unstable directions contribute to diffusion only when activation is not used, that is when 
mechanism B (rather than A) rules the dynamics. Therefore, one may conclude that the 
approach of p2] works so well exactly because above Tx activated processes are not relevant. 
It could be potentially very useful to reconsider the role of INM unstable modes in the light 
of our scenario ||23l . 
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